Exponential super-twisting control for nonlinear systems with unknown polynomial perturbations

The study focuses on the control of nonlinear dynamic systems in the presence of parameter uncertainties, unmodeled dynamics, and external disturbances. The lumped perturbation is assumed to be bounded within a polynomial in the system state with the polynomial parameters and degrees unknown a priori such that it accommodates a quite wider range dynamic systems. Based on the studies in recent super-twisting algorithm designs and the idea from adaptive sliding mode control for nonlinear systems with uncertainties, we propose a novel adaptive super-twisting algorithm with exponential reaching law, or exponential super-twisting algorithm (ESTA), for the high-stability and acceptable accuracy control of the aimed nonlinear dynamics. The stability analysis and practical finite-time (PFT) convergence are proven using Lyapunov theory and an intuitive analysis of the control behaviour. Simulations are performed to compare the proposed ESTA with the existing super-twisting method and the traditional proportional integral differential control. The simulation results demonstrate the effectiveness of the proposed ESTA in terms of the fastest settling time and the smallest overshoot.

To overcome the drawback of the overestimation and to deal with uncertainties of unknown bounds, the design of adaptive sliding mode control (ASMC) was introduced 2,[34][35][36][37][38] where the time-varying switching gain is designed to adaptively compensate for the lumped uncertainties.Most ASMC techniques use integral adaptation laws or integral adaptation law combined with other techniques such as the σ-modification 34 or the dead zone method 35 .It is shown that, by using integral-type adaptation law, the systems with uncertainties of unknown constant bounds have finite-time stability 36 .However, the system response to the perturbation is relatively slow and, even though the overestimation is avoided somehow, the chattering phenomenon is still observed 39 .To further attenuate the chattering phenomenon, a possible adaptation law is to reduce the switching gain to a minimum admissible value since the magnitude of the chattering level is proportional to the magnitude of the switching gain 37 .In 37 , the ASMC law is applied to STA and uses a low-pass filter to tune the switching gain in the control to a possibly minimum value once the sliding mode is established.However, it still requires a large enough feedback gain to compensate the perturbation with affine function form in the system state, and the use of a low-pass filter introduces a time delay which affects the transient phenomenon.
To achieve fast response and chattering-attenuation property, asymptotic reaching laws rather than ASMC can be found as the form of power reaching law 1 and an exponential reaching law 40,41 .However, they either lose the robustness when system states are around the sliding surface 1 or require a priori knowledge of the uncertainty bounds 40,41 .Yang et al. 42 proposed an ASMC technique using exponential reaching law, called adaptive exponential sliding mode control (AESMC), for the control of a beaingless induction motor.Still, for the AESMC design in 42 all the bounds of the uncertainties must be known a priori.In brief, the aforementioned STA techniques assume that the bounds of the uncertainties are known a priori, while the ASMC techniques assume that the uncertainties are bounded within some constants or affine functions in the system state with the bounds a priori known or unknown.For the case that the bounds of the uncertainties are unknown a priori, seldom researches studied the STA based ASMC.More than the case that the uncertainties are bounded within constants or affine functions in the norm of the system state, the uncertainties may be bounded within some quadratic or cubic functions in the norm of the system state 3 with the bounds unknown or inaccurate a priori.For instance, with inaccurate parameters the aircraft dynamics and robot manipulators contain uncertainties bounded by quadratic functions of the state.Also, polynomial uncertainties can be found in the systems of the inverted pendulum mounted on a cart, the tunnel-diode circuit dynamics and the Duffing dynamics 3 .
Such cases (with quadratic or cubic uncertainties) represent a wider class of nonlinear systems than the above STA and ASMC techniques considered.More generally, all the above discussed uncertainties can be assumed to be bounded within some polynomials in the norm of the system state with the bounds unknown a priori.To stabilize the above wide class of nonlinear systems, we propose a control method of an exponential super-twisting algorithm (ESTA) where the algorithm structures from STA are integrated with a novel exponential reaching law.Thus, the main contributions of the work can be illustrated as follows.
• Compared to the conventional STA where the bounds of the system's uncertainties must be known a priori, or the adaptive STA where the uncertainties are assumed to be bounded within some constants or affine functions, a more general case of nonlinear systems is considered in this study where the uncertainties are assumed to be bounded within some polynomials in the norm of the system state with the bounds unknown a priori, i.e., both the polynomial parameters and degrees are unknown a priori.• A novel algorithm, adaptive super-twisting algorithm with exponential reaching law, is proposed to stabilize the aimed nonlinear systems.The stability and the practical finite time convergence 43 of the new design are proven using Lyapunov theory and an intuitive analysis of the control behavior.• The proposed new ESTA is compared with the traditional STA and PI methods in simulations, and the simulation results demonstrate the superiority of the new design.
Section "Problem formulation" introduces the nonlinear systems with unknown polynomial uncertainties.The STA technique is recalled with the discussion of the stability issues in some STA designs over past decade.The new ESTA laws are proposed in section "Methodology".The new design with single-input-single-out (SISO) form is introduced first for the ideal case.Then the multi-input-multi-output (MIMO) forms of the new ESTA are designed for the ideal and real cases.Simulation results are presented in section "Simulation results" to verify the effectiveness of the proposed algorithms.Finally, section "Conclusion" concludes the paper.

Problem formulation
In this section, we recall the existing STA techniques 9,10,16,26 to analyze the stability limitations.Useful definition and mathematical lemma are first introduced in this section.

Definition 1
The signum function is given as Then, Lemma 1 If the time derivative of σ (t) exists, then

Problem statement
Consider the uncertain nonlinear dynamic system where σ ∈ χ ⊂ R n is the measured signal designating the system state or any sliding variable, t ∈ R + is the time, and u ∈ R n is the control input signal.
0 is an equilibrium of (4).Function f (σ , t) ⊂ R n represents lumped perturbation containing parameter uncertainties, unmodeled dynamics, and external disturbances and Function g(σ , t) ⊂ R n×n contains parameter uncertainties.
Assumption 1 The norm of the perturbation f (σ , t) is upper-bounded with some unknown polynomials in the state vector σ ∈ χ .More specifically, where q is an uncertain finite integer, a i ( i = 0, 1, • • • , q ) uncertain non-negative finite values, and n i positive real scalars.
Note, Assumption 1 includes but not limits to the following example.

Example 1
Assumption 2 Let g(σ , t) T be the transpose matrix of g(σ , t) .The uncertain term g(σ , t) is positive definite in wide sense, i.e., its symmetric part g s defined by is positive definite in the regular meaning.
The assumption 2 implies that the minimum eigenvalue of g s is lower-bounded by a positive finite constant.In other words, there exists a positive finite constant b such that where I n is an identity matrix with dimension n.
Remark 2 Usually the control parameter g(σ , t) is assumed to be known a priori.For the case where g(σ , t) is varied, the bounds of the variations usually are assumed known a priori as well.In this contribution, we consider the statement (7), that is, the variation values of the control parameter g(σ , t) are loosely structured, which is another way extending the class of nonlinear systems to be addressed.
(3) www.nature.com/scientificreports/ The selections of the scalars k 1 and k 2 in the existing STA designs are based on some stability analyses and have different complicated forms 9,16,26 .Alternatively, k 1 and k 2 are simply selected 10,11,15,18,19 as In 10,11,15,16,18,19 it is stated that, if the scalars k 1 and k 2 are sufficiently large or respectively satisfy ( 12) and ( 13), the finite-time stability of the STA design ( 9) is guaranteed.
The condition (11) contains a possible strict ramp disturbance in the domain of time t.Consequently, the magnitude of the ramp disturbance may go to infinity as time elapses.However, most the existing STA control designs are difficult to handle a disturbance with extremely large (or infinity) amplitude.It is possible that some researchers are aware of the inappropriateness or unreality of this hypothesis.The condition (11) is not presented in STA designs of 20,22,45 .Papageorgioua and Edwards 19 used a restricted assumption, a T-periodic perturbation with (11), to replace the simple (11).Base on the 'T-periodic' assumption, the authors further demonstrate that under smaller gain conditions ( k 2 < L ), the solutions of the closed loop system converge to a stable limit cycle around the origin as well.
Remark 3 Note, the Assumption 1 includes the condition (10) but not (11).In real dynamic systems, the condition Eq. ( 11) which possibly contains a strict ramp disturbance is unnecessary and has little practical significance.
Based on the above discussion, we take Assumptions1 and 2 instead of the conditions (10) and (11) in the following new ESTA design.

Methodology
In this section, we will propose a novel adaptive STA based control law, ESTA design, to expand the traditional STA design by using an exponential reaching law to constrain the system states of ( 4)- (7) to the vicinity of zero in finite time.The finite time convergence to the vicinity of the equilibrium is also defined as practical finite time stability (PFTS) in 43 or finite settling time stability (FSTS) in 46 .We start the new ESTA design for the system (4)- (7) for the ideal scalar case.Then we extend the new design to the real case with MIMO form.

Remark 4
The PFTS means that the state σ in (4) converges to the vicinity of the equilibrium in finite time.The finite settling time stability (FSTS) defined in 46 has a similar meaning.

ESTA design for SISO systems
We first consider an simple SISO case where the dynamic systems are ideal 9 .Consider the following control law for some positive constants c 1 and c 2 .
Proof Consider the situation where the system state σ is outside of the vicinity of the equilibrium, i.e., for any small positive value of ε .We will prove that the state will converge to the domain |σ | ≤ ε in finite time.Note from (5)   Since the exponential term bc 1 e √ |σ | will eventually be bigger than the polynomial term bc 16) is upper bounded, i.e., (10) Vol.:(0123456789) Now we consider two cases of σ .Case I, the system state is in the positive space of the vicinity of the equilibrium, i.e., σ ≥ ε > 0 .We obtain from ( 17) and ( 18), Noting that c 2 b(t − t * ) is continuously increasing as time elapses after t ≥ t * .That is, the rate of the decline of the state |σ | is getting faster and faster as time elapses and the state |σ | will eventually reach the region |σ | ≤ ε in finite time.By integrating (20) between t * and t > t * and using the Comparison Lemma 3 , we obtain Then the reaching time t F is estimated as Case II, the system state is in the negative space of the vicinity of the equilibrium, i.e., σ ≤ −ε < 0 .We obtain from ( 17) and ( 18), Note that the inequality ( 22) has a same form as (20) in Case I. We conclude that |σ | will eventually reach the region |σ | ≤ ε in finite time and remains on it thereafter with the reaching time estimated as Theorem 1 shows that for any positive constants c 1 and c 2 the control (14) will always steer the system state to the vicinity of the equilibrium point in finite time even if the lumped perturbation contains some polynomial form in the system state and even if we do not known the boundaries of the uncertainties a priori.The reaching time or the finite settling time t F of Theorem 1 contains two parts, the compensating time t * and the converging . One can see that the settling time t F depends on the control gain c 2 and system parameter b .Generally, t F can be reduced by increasing c 2 and b values.An insight discussion of the finite settling time t F estimation can be found in Ref. 39 .
( www.nature.com/scientificreports/One drawback of the conventional sliding mode control is the chattering phenomenon where the switching function sgn(•) is the main source of it.In real implementation, the function is usually replaced by some smooth functions to attenuate the chattering effects 2,3,37 .In this work, a simple smooth function sgn(σ ) ≈ σ |σ | + µ with a small positive scalar µ 47 can be used to replace the switching function.
where the small positive constant µ is related to the thickness of the boundary layer of the real sliding mode 3 .
Using the aforementioned small positive scalar ε as the vicinity of the equilibrium, we have the following proposition.
The proof is similar to the proof of Theorem 1 and is shown in Appendix "Proof of the Proposition 1".

Remark 5
The smoothing function can be used in most sliding mode control.The main imperfection of the smoothing function is that it is not suitable for direct application of a small number of SMC controls, such as the STA design (9).In other words, when the smoothing function sgn(σ is chosen to approximate the switching function, the existing STA design (9) with conditions ( 10) and ( 11) may encounter stability issues which can be seen in the following example.
The above example shows that by using the existing STA design, if the switching function is approximated by some smoothing functions, the system may become unstable no matter the gain k 2 satisfies (13) or k 2 is sufficient large.(See Fig. 1 for k 2 = 1.1 satisfying (13) and Fig. 2 for k 2 = 10 sufficiently large.)

ESTA for MIMO systems
Now we extend the new ESTA design for the MIMO systems (4)- (7).Let σ be the 2-norm of the system state σ .Consider the following control law for some positive constants c 1 and c 2 .
Theorem 2 Consider the MIMO system (4)-( 7) subject to (24), then the state in (4) has PFTS.www.nature.com/scientificreports/Proof In the following, the argument t of most given vector field (i.e., σ , f, g, etc.) will be omitted for simplicity.Consider the following Lyapunov function candidate Let u T be the transpose vector of u in (24).Using ( 4) and ( 24), the time derivative of V along the system trajec- tories is Note, the scaler Then, we obtain from ( 7), ( 26) and ( 27), Consider the situation where the system state σ is outside of the vicinity of the equilibrium, i.e., �σ � ≥ ε .Using �σ � = √ V , �σ � 2 = σ T σ , and recalling Assumption 2, we reformulate (28) as, (25) The inequality (29) can be rewritten as For any �σ � ≥ ε , the term c 2 b t 0 σ (τ ) dτ in (31) keeps increasing and will eventually compensates for h .Since this compensating action occurs for any �σ � ≥ ε , we conclude that there exists a time instant t * such that h ≤ c 2 b t * 0 �σ (τ )�dτ .Then, after the time instant t * , i.e., t ≥ t * , We obtain from (31)   By applying the mojorant curve approach 39 , we conclude that σ reaches the region �σ � ≤ ε in finite time and remains on it thereafter with the reaching time estimated as The above ESTA algorithm ( 24) is designed for the ideal case.For the real implementation, the magnitude of the integral term t 0 σ (τ ) dτ may rise to a very large value since it keeps growing for any �σ � � = 0 .Also, the term σ as a denominator in (24) may cause some singularity problem when σ is close to zero.To eliminate these problems, we use ε-adaptation 36 and a smoothing function 47 in the following ESTA design for the MIMO systems ( 4)-( 7) in real implementation.
where c 1 and c 2 are positive gains and µ is the smoothing factor in the smoothing function σ �σ � + µ used for replacing the switching function sgn(•).
The proof is ignored here to save the space.
Remark 6 Note that the controllers ESTA (14) and STA (9) use full state feedback control.These controllers can be applied to higher-order systems if they can be converted to first-order systems and the state of the firstorder system is available.If the state of the system cannot be measured directly, but the system is observable, it is often possible to construct an "observer" or simply use some differentiator to estimate the full state.When the "observer" or differentiator are good to use, the two-part structure of proposed ESTA methods ( 14), ( 23), ( 24) and ( 33) can be extended to three-part structure.For example, ESTA controller ( 14) can be extended to a controller containing a differential term.
where |σ | represents the estimation of the differential of the measured sensor signal σ if σ is not directly available.The stability proof of the control (34) is verbose and ignored here.
The selection of gains can refer to the gain adjustments of PID control.Generally, the exponential gain c 1 in proposed ESTA methods can be first tuned to a relatively large value.Then gain c 2 can be tuned to eliminate the steady error.If the measured signals are differentially reliable, c 3 can also be tuned to further reduce the overshoot.

Remark 7
The reference signal of the controllers designed in this paper can be time-invariant, such as coordinate values, or time-varying, such as trajectories.For the trajectory reference signal, we require it to be second-order differentiable such that the error dynamics equation of the trajectory tracking control system can be constructed.

Simulation results
Although the STA algorithm appeared in the 90s of the last century, it is still widely used today, such as the literatures 13,14,19,25,[29][30][31][32] of the last three years.The core structure of the STA algorithms used in these literatures is consistent.As a result, these newer STA control methods also share common features, such as a reduced settling

Illustrated example
In this section, an example which contains external perturbations and polynomial uncertainties in the norm of the state is given to illustrate the system response of the proposed ESTA.Consider the system (4) with g(t) = b = 1 for simplicity.The lumped perturbation f (σ , t) = f 1 (t) + f 2 (σ ) is chosen as where f 1 (t) represents external disturbances and f 2 (σ ) contains system parameter uncertainties or unmodelled dynamics.In particular, α 3 (σ 2 + σ 3 ) represents higher order polynomial bounded disturbances if the unknowm scalar α3 = 0 .In simulation the proposed ESTA control (23) (for real case) is compared to the existing STA control ( 9) and the conventional proportional-integral (PI) control which has the following form: Note, the three methods have similar proportional-integral forms.For simplicity and comparability, the 'P-parameters' in ( 9), ( 23) and ( 36) are selected as k 1 = c 1 = d 1 = 1 , as well as the 'I-parameters' k 2 = c 2 = d 2 = 1 .We choose the smoothing factor µ = 0.001 and keep it unchanged for all simulations.

Transient response under external disturbances
To test the system transient response under various external disturbances, we let the parameter of the unmodelled dynamics α 2 = α 3 = 0 be fixed.Two levels of the external disturbance, α 1 = 2 and α 1 = 5 , are used to respectively represent a moderate and a relatively large external disturbances.We first choose a moderate external disturbance with α 1 = 2 where the simulation results of the state σ (t) and the control input u(t) are shown in Figs. 3 and 4, respectively.Then we increase the magnitude of α 1 to a relatively large value, i.e., α 1 = 5 where the simulation results are shown in Figs. 5 and 6.One can see that, compared to the STA and PI methods, the proposed ESTA method has the fastest settling time and the smallest overshoot.In particular, for the moderate external disturbance with α 1 = 2 , the system settling time t s of the proposed ESTA is reduced from 10 and 4 seconds to 3.2 second ( 68% reduction and 20% reduction) compared to those of the existing PI and STA methods, respectively.For the relatively large external disturbance α 1 = 5 , the t s of the pro- posed ESTA is reduced from 8.1 and 9.2 seconds to 6.1 second ( 25% reduction and 34% reduction) respectively compared to those of the PI and STA methods (see t s in Table 1).Simultaneously, compared to the PI and STA methods the percentage of the system overshoot O.S. by using the proposed ESTA is respectively reduced from 30% and 40% to 25% for the moderate external disturbance.Meanwhile, the O.S. by ESTA is reduced to 16% for the relatively large external disturbance (see O.S. in Table 1).Note, the percentage of O.S. is calculated from the control signals (Figs. 3 and 5) because we choose the system steady state at its equilibrium state, i.e., σ (∞) = 0.

System response under external disturbances and linear unmodelled dynamics
To test the system response under various linear parameter uncertainties (or unmodelled dynamics), we keep the magnitude of the external disturbance unchanged, i.e., α 1 = 2 , no higher order polynomial disturbance, i.e., α 3 = 0 , and choose different values of linear unmodelled dynamics α 2 .
First we choose a relatively small value of α 2 = 0.5 .One can see that the ESTA method still has the fastest settling time and the smallest overshoot (see Figs. 7 and 8).Specifically, the settling time t s dropped from 18 and 5.2 seconds respectively by applying PI and STA approaches to 3.2 seconds by using the proposed ESTA method (see t s in Table 2).Then we increase the parameter uncertainty to a moderate level of α 2 = 1.0 .From the simulation results Figs. 9 and 10, one can see that by using STA and PI methods the systems are unstable and oscillating, respectively.In contrast, the system is still stable by using the proposed ESTA control method.
Note that any α 2 > d 1 = 1 will lead to an unstable PI control as well.We continue to increase the α 2 value to 1.5 where the simulation results are shown in Figs.11 and 12.It can be seen that the system becomes unstable by   36) (dash red), ESTA ( 23) (dot blue), and STA (9) (solid black).
Table 1.Settling times (obtained from Figs. 4 and 6) and overshoots (obtained from Figs.    www.nature.com/scientificreports/using the PI method as well as by using the STA method.In contrast, the system state is still stable by applying the proposed ESTA method.

System response under high order polynomial bounded disturbances
Now we test the system response where higher order polynomial bounded disturbances arise.We keep the magnitude of the external disturbance and linear unmodelled dynamics unchanged, i.e., α 1 = 2 and α 2 = 0.5 .
We choose different values of α 3 to see the system response when higher order polynomial disturbances arise.First we choose a relatively small value of α 3 = 0.1 .One can see that the ESTA method still has the fastest set- tling time and the smallest overshoot (see Figs. 13 and 14).Specifically, compared with PI and STA approaches, the settling time t s dropped from 15 and 3.1 seconds to 2.3 seconds and the overshoot O.S. reduced from 90% and 70% to 50% by applying using the proposed ESTA method, respectively (see t s and O.S. in Table 3).Then we increase the magnitude of polynomial disturbance to a moderate level of α 3 = 0.25 .From the simulation results Figs. 15 and 16, one can see that by using the PI method the systems become unstable.When the magnitude of the polynomial disturbance coefficient is increased to α 3 = 0.57 , the STA method starts to become unstable as well (see Figs. 17 and 18).In contrast, the system keeps stable by using the proposed ESTA control method for the  above three cases.Moreover, the system is still stable for ESTA as the coefficient value of polynomial disturbance arises to α 3 = 1 (see Figs. 19 and 20).

Inverted pendulum
An inverted pendulum mounted to a motorized cart is commonly founded in control system textbooks and research literature 48 .In Fig. 21, M is the cart mass, m the pendulum mass, 2l the length of the pendulum, and b the coefficient of friction of the cart.In this case we only consider a two-dimensional problem where the pendulum is constrained to rotate in the vertical plane.The control input is the force F that moves the cart horizontally and the output is the angular position θ for simplicity.Since the angular velocity of the pendulum can be estimated via a filtered differentiator of the angular position θ , a PID rather than a PI controller is often used in most literature to stabilize the pendulum system 48 .For the comparison, we stabilize the system using ESTA, PID and STA methods.

Nonlinear equations of motion
Applying Newton or energy method, we can derive the following nonlinear equations of motion.) by using three control methods for the nonlinear system (4) under higher order polynomial disturbances: α 1 = 2 fixed, α 2 = 0.5 fixed, and α 3 = 0 varied.Significant values are in [bold].where I = (1/12)m(2l) 2 = 1 3 ml 2 is the mass moment of inertia of the pendulum around the center of the pen- dulum and g is the gravity.Then we can solve the above equations ( 37) and (38) to obtain the following full nonlinear equation.  36) dash red), ESTA (23) (dot blue), and STA (9) (solid black).One can see that the PI method becomes unstable.www.nature.com/scientificreports/where q = 4 3 l(M + m) − ml cos 2 θ > 0 .Let σ 1 = kθ + θ be a sliding variable 10 and u = −F be the new control input.The second-order differential equation (39) can then be converted into a first-order differential equation.

Stabilization of the inverted pendulum
Now we can build the simulation model of the inverse pendulum system using the equations of motion (39)  and (40) and applying the three control laws (34), (43) and (44) for the three methods ESTA, PID, and STA, respectively.The parameter values for the inverted pendulum system are shown in Table 4.
We test the transient response of the inverted pendulum system under different pulse disturbances.First we tuned the PID control, and, after some try and error, it is found that the PID gains [k 1 , k 2 , k 3 ] = [50, 10, 20]  provides a satisfactory response.For comparison, let's make the corresponding gains the same for all three methods, i.e., k      , the overshoot and settling time of the system are significantly different between the three methods (Figs. 23 and 24).In particular, for the pulse angle offset A o of medium values 22.9 • and 45.8 • , the system settling time t s of the proposed ESTA is reduced from 3 and 0.8 seconds to 0.4 second ( 83% reduction and 50% reduction) and 5 and 1.2 seconds to 0.9 second ( 81% reduction and 25% reduction) compared to those of the existing PI and STA methods, respectively (see t s and O.S in Table 5).For A o of large values 68.8 • and 80.3 • , one can see that by using STA and PI methods the systems are unstable.In contrast, the system is still stable by using the proposed ESTA control method (Fig. 24).Note that an offset angle A o ≥ π 2 = 1.57radius ( 90 • ) means that the inverted pendulum would fall down com- ) can be a possibly stabilizable range of the inverted pendulum.The perfor- mance of the inverted pendulum stabilized via the three methods is summarized in Table 5.

Phase portraits of angular positions and angular velocities
The simulation results of the inverted pendulum system also let us analyse the phase portraits of the system.In the phase portraits Figs.25-27 the angular position state θ versus the angular velocity state θ are plotted along with two levels of local magnifications.For the system under pulse disturbances with relatively small offset angles ( A o = 0.1 and A o = 0.2 , Fig. 25) and medium offset angles ( A o = 0.4 and A o = 0.8 , Fig. 26), one can see that the trajectories of the two phases θ and θ have similar performance with the three methods.All θ − θ trajectories start vertically upwards from the origin and eventually return in a spiral towards the origin.Compared to the STA and PID methods, the trajectory ranges of the ESTA method are slightly larger, but the returns are much closer to the origin.For the offset angles increased to large values ( A o = 1.2 and A o = 1.4 ), the proposed ESTA method performed much better than the STA and PI methods (Fig. 27).Specifically, not only the θ − θ trajectory regression to the origin of the ESTA method is faster and closer, but also the trajectory range is much smaller.

Conclusion
The paper first investigates some perturbation assumptions in some STA designs over the past decade.The polynomial perturbation is then used to fix the stability issue and further extend the range of nonlinear dynamic systems.To handle the extended systems, a novel adaptive super-twisting sliding mode control with exponential reaching law is proposed.The new ESTA method is extended from the existing STA design and integrated with an novel exponential reaching law.The stability analysis and finite-time convergence are proven using Lyapunov theory and an intuitive analysis of the control behavior.The new design successfully applied to the control of the nonlinear systems having lumped perturbations bounded by a parameter-unknown-a-priori polynomial.The simulation results of the illustrated example show that, for the external disturbance, the system settling time t s of the proposed ESTA is reduced by 25 − 68% and 20 − 35% compared to those of the existing PID and STA methods, respectively.Simultaneously, compared to the PI and STA methods the percentage of the system overshoot O.S. by using the proposed ESTA is respectively reduced from 30% and 40% to 16 − 25% .It is worth noting that when the polynomial disturbance appears and arises to medium or relatively large amplitudes, the system using STA and PI methods began to become unstable, but with the proposed ESTA control method it remains stable.The simulation results of the inverted pendulum also show that, when the inverted pendulum  is disturbed by small pulses close to the equilibrium point, the use of all three methods has a good stabilizing effect.However, when the pulse is of medium size, ESTA has better stability than PID and STA.Specifically, ESTA method reduced settling time by 25 − 80% and overshoot by 10 − 50% .Moreover, when pulse is large, ESTA can still stabilize the system quickly, while PID and STA cannot stabilize the system.Compared with the commonly used STA and PID methods, the proposed ESTA approaches exhibit superiority in terms of faster settling time, smaller overshoot and stronger stability.Moreover, ESTA methods are able to deal with nonlinear systems where higher order polynomial bounded disturbances may arise.Thus, the proposed ESTA methods are of great significance for future designs and applications of robust control for a large class of nonlinear dynamic systems, such as all-weather flight control of civil airliners, flight control of rescue helicopters in disaster weather, motion control of humanoid robots subjected to sudden external forces, etc.. Our future work will focus on a real flight control of an unmanned aerial vehicle in strong wind environments.

Appendicies Proof of the Lemma
Proof There are two cases of |σ | � = 0.
Note, the above lemma does not include the case of σ = 0 where the system is under the steady state.

Figure 1 .
Figure 1.Simulation results of the Example (2) by using the existing STA control (9) with the switching function sgn(σ ) replaced by σ |σ | + 0.001 for the case I k 2 = 1.1.

Figure 2 .
Figure 2. Simulation results of the Example (2) by using the existing STA control (9) with the switching function sgn(σ ) replaced by σ |σ | + 0.001 for the case II k 2 = 10.

Figure 21 .
Figure 21.Schematic representation of the inverted pendulum system.

Proposition 1 ProofFigure 27 .
Figure 27.Phase portraits of angular positions and velocities under pulse disturbances with relatively large offset angles 1.2 radius ( 68.8 • ) and 1.4 radius ( 80.3 • ) using ESTA (23) (dot blue), PID (36) (dash red), and STA (9) (solid black). /b 1 sgn( |σ |) , but a decrease in robustness.Compared to the STA methods, the advantage of the ESTA proposed in this paper is that it maintains a short or shorter settling time, while at the same time has a significant improvement in robustness.Because of the comparison in terms of settling time (stability time) and robustness, we chose two quantitative indicators, settling time t s and overshoot O.S., for stability analysis.

Table 4 .
Parameter values for the inverted pendulum mounted to a cart.

Table 5 .
Settling times and overshoots obtained from Figs. 22, 23 and 24 by using three control methods for the inverted pendulum system, (39) and (40), under external pulse disturbances with pulse width 0.2 seconds.Significant values are in[bold].